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A COMPUTER SIMULATION STUDY OF MARK-RECAPTURE 
METHODS IN ECOLOGY 


1. Means and empirical standard deviation from four models 


Ву С. ETTERSHANK* AND DAPHNE L. ETTERSHANK* 


ABSTRACT: А Fortran computer program is described which simulates an animal 
population (with births and deaths) subjected to mark-recapture study. All parameters 
of the study can be specified. The population is estimated by the models of Lincoln, 
Schnabel, Bailey and Jolly, and mean values and empirical standard deviations are 
derived on the basis of twenty replications of each set of experimental parameters. 


INTRODUCTION 


The use of marks on animals to aid the study 
of their biology is a long-established one. Izaak 
Walton, in 1653, relates that Sir Francis Bacon 
attached ‘ribands’ to young salmon so they could 
be identified on their return. Marking is often 
used for the estimation of population size, and a 
great deal of sophistication has now been brought 
to the mathematical models used in making such 
estimates (Cormack 1968). 

Considcrably less attention has becn given to 
the testing of these models on known populations: 
such studies have usually involved the comparison 
of several models on a single population, as, for 
example, Parr's (1965) study of a population of 
dragonflies using five models. Phillips and Camp- 
bell (1970) used four models on a closed popula- 
tion of whelks on a rocky shore platform. 

Manly (1970) simulated populations subjected 
to mark-recapture methods on a computer, and 
obtained data on the behaviour of four estima- 
tion modcls. Later (1971a) hc examined the 
bchaviour of the varianee estimates produced 
by the method of Jolly (1965) by simulation 
methods. More recently, other aspccts of Jolly's 
method have been examined (Cormack 1972, 
Manly 1971b, Carothcrs, in prcss), and these 
are discussed more fully below. A very use- 
ful bibliography of mark-recapture methods has 
been compiled by Anderson (1972); compre- 
hensive discussions are given in Robson (1969) 
and Scber (1972). White (1971a & b) supplies 
a Fortran computer program which will handlc 
estimation by Jolly's (1965) method. 

Models of mark-recapture estimation are bascd 


on a series of assumptions about the population 
under study. The accuracy of the estimate will 
depend to a large extent on how closely these 
assumptions are met, and on the sensitivity of 
the particular model to each individual assump- 
tion. Thus apart from sampling errors, thcre are 
also biological errors in the estimation process. 
Thc assumptions are as follow: 


1. The marking method does not impair the 
animal functionally or behaviourally. 

2. Marks are permanent, unalterable, and accu- 
ratcly decipherable. 

3. The return of marked animals into the popula- 
tion does not alter its spatial structure. 

4. АП animals are equally available for capture, 
irrespective of 
(a) mark status (i.e. whether marked or un- 

markcd) 

age, sex or other group existing in the 

population (whether distinguishable ог 

cryptic) (ie. each group is sampled in 
the same proportion as it occurs in the 
population) 

(c) previous eapture history (i.e. being сар- 
tured neither increases nor decreases the 
probability of an animal's being recap- 
tured subsequently). 

5. All parts of thc habitat have an equal chance 
of being sampled. 

6. Sampling is carried out in a diserete period 
of time which is small compared with the 
total time of the study. 

The following assumption is also said to hold 
for the Lincoln Index and Schnabel models: 

7. The population is either closed, i.e. there is 


(b) 
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neither recruitment (births and immigration), 
nor losses (deaths and emigration) to the 
population of the period of study; or clse 
recruitment and/or losses can be independently 
assessed and allowed for. 


Assumptions 1 and 2 can be examined by inde- 
pendent laboratory and field studies. It does seem 
probable, however, that some marking methods 
must impair the recipient: for example the more 
radical toe-removal methods uscd to code small 
mammals. Such methods are undoubtedly perma- 
nent, and unless the animal accidentally loses 
another toe, are accurately decipherable. Paint- 
marks used on inseets can chip off, or colours 
can fade. Such marks ean be field-tested by using 
two (or more) separate, independent marks on 
each animal, thus checking on the long-term relia- 
bility of each system. Marking systems should 
not be of such complexity that human error in 
reading marks becomes a factor. This can be 
assessed by having more than one observer inde- 
pendently decipher each mark. 

Assumption 3 is a procedural matter, in the 
hands of the researcher—it requires deep thought 
into the timing, as well as the placement, of 
releases. It is probably also advisable to release 
animals at the point of capture, and not to rely 
on dispersal from one or a few release points. 

Before examining the other assumptions more 
fully, it is necessary to consider the process of 
estimation (Cormack 1972). The population con- 
tains a sub-population of identifiable individuals, 
the members of which bear either unique, indi- 
vidual marks or a code-mark which identifies the 
date (and perhaps other information such as 
place) of marking. It is assumed that the charac- 
teristics of this sub-population are known to the 
biologist (e.g. its size), and that any changcs are 
also known or can be detected (e.g. additions and 
deletions from the sub-population). Most import- 
antly, it is assumed that changes estimated from 
this sub-population of identifiable individuals arc 
an estimate (whose accuracy can, in turn, be 
estimated) of processes in the population. Thus, 
for each capture-reeapture cycle, estimates are 


£ р‹ = т/М, 


made of sampling intensity р; (in the notation of 
Cormack 1968)*, birth rate B,, and probability of 


survival фі. 
From these, further estimates may be made 
for the overall population, e.g. №, = nM m; 


(= n,/p,)*. Cormack (1972) points out that to 
solve the estimation formula it is necessary to 
know (or to directly estimate) one and only one 
of the sets Мь р, №, Ф; or В;; knowledge of the 
set Ф or the set В; must be supplemented by 
knowledge of one М; or N,. From this logical 
base all the models of mark-recapture estimation 
can be interpreted. Thus Jolly's (1965) and Seber's 
(1965) t maximum likelihood solution of the esti- 
mation of the marked sub-population is 
М, = Ju -- 5; 2:۳ 

where 2; is the number of animals marked before 
the time of the ith samplc, but not included іп it; 
and r; the number of the s; animals released at 
the ith period that are scen later. Thus the second 
term is an estimate of the number of marked 
animals alive but unseen at the ith sampling 
period, and т, is the number аешаПу seen (Jolly 
1965). 

On the question of availability of animals for 
capture, detailed observation may show whether 
this is likely to be violated. If certain individuals 
or groups are more conspicuous in the population, 
their probability of capture may be biased (posi- 
tively or negatively). Conspicuousness may be 
physical (e.g. bright plumage or bizarre ornamen- 
tation on birds or insects, as in sexual dimor- 
phism), or bchavioural (e.g. due to dominance 
hierarchies or division of labour (Golley & Gentry 
1964, Wilson 1971) ). Cormack (1972) shows that 
the catchability of these individuals is higher at 
the start and thus average catchability declines 
during the course of the project. However, if the 
catching method discriminates against conspicuous 
individuals (e.g. by personal bias of an operator, 
ог Бу a size-filter effect of a trap), the reverse or 
a different bias may arise. 

Another aspect of this same question is trap- 
ping history, and trap-proneness and trap-shyness 


та is the number of individuals of the sub-popnlation of marked animals М, sampled оп the ith 


where 
оссазїоп 
$&‹=М, + 1/¢ М.+я— т.) 
where 5, is the number of marked animals released after the ith sample. 
B, = А, + 1-- (К. + و‎ — па) 
where л; is total sample size оп the ith occasion, № is the estimated population at time i. 


1 Jolly (1965) and Seber (1965) produced the same solutions to the problem of mark-recapture estimates of 
populations, under a stochastic model; Jolly’s treatment, however, is more generalized, and is referred to 


hereafter. 
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are well known to workcrs on vertebrate animals. 
(The effects of training on invertebratcs should 
not be overlooked, howcver.) Orians and Leslie 
(1958) provide a test for random recapture of 
individually-marked animals. One further situation 
to be avoided is the use of baited traps in, for 
example, studies of social insects where rccruit- 
ment of a specialized group of foragers occurs 
when а single suceessful forager has returned to 
the nest (Wilson 1971). The effects of bias in 
catchability is reflected in the estimate of М; this 
is discussed for the Jolly model by Cormack 
(1972) and Carothers (in preparation). 

The bias in estimates caused by differences in 
survival is discussed by Cormack (1972). In the 
Jolly method, general survival is estimated for 
each period; in others (e.g. Fisher & Ford 1947) 
an Overall estimate is made for the entire pcriod, 
and this may lead to serious bias. Thc most serious 
effect of differential survivals will occur when 
this effect is associated with marking. The situa- 
tion can be allowed for if the effect is assumed to 
occur only after the first capture (and marking), 
and to last only until the next recapture (the type 
I loss of fisherics biologists) by considering two 
sub-populations of newly-marked and mature- 
marked individuals. The cstimates then follow 
logically (Cormack 1972). Manly (19715) dis- 
cusses the same situation for the case where ani- 
mals bear individual marks—he develops a method 
for analysing type I losses, as a spccial case of 
Jolly’s cstimator (if the animals die while being 
marked, this is allowed for in the Jolly model). 

Assumption 5, that all parts of the habitat have 
an equal chance of being sampled, rests on the 
design of the sampling program, which is excel- 
Ісі йу covered in specialized texts. Planning for 
efficient sampling will also ensure that assumption 
6 is met—that sampling is carried out in a rela- 
tively short, discrete time intcrval. 

The present paper presents the results of a 
computer simulation study in which four cstima- 
tion models were examined. The sclection of these 
models was influcnced by thcir being used or 
under consideration for use in studies of the Desert 
Biome Program of the US/IBP Analysis of Eco- 
systems. As mark-recapturc studies were to be 
used for a wide range of taxa, a comprehensive, 
realistic model of a population subject to such 
study was required. This paper describes the com- 
puter model and the first sct of results obtained 
from it, for idealized populations which com- 
pletely meet the assumptions discussed above. 


METHODS 


Idealized animal populations, being studied by 
mark-recapture methods, were simulated by a 


program written in Fortran IV called SYNPOP. 
This was run initially on an IBM360/50 eomputer 
at New Mexico State University, Las Cruces, New 
Mexico, and later adapted and expanded for a 
Burroughs B5500 computer at Monash University. 

The basic outline of the program is relatively 
simple (sce Fig. 16). An array ‘POP’ is set up, and 
a number of places in the array—corresponding 
to the starting or ‘seed’ value of the population— 
are set equal to 99, which identifies new-born 
animals. They retain this value unless they are 
captured (see later) or killed (when their value 
becomcs zero). 

The population then goes through a series of 
cycles of dcaths, births and sampling, the number 
of cycles being specified by the experimenter. He 
also spccifies death and birth rates, either as 
probabilities or as fixed fractions of the population 
at the start of each day, and these may apply 
uniformly throughout the period of the experi- 
ment, or may be specified (independently) in 
particular cycles. The experimenter may also spe- 
cify that births should exceed deaths or vice versa; 
and finally, he may nominate a probability that 
animals captured in samples, are killed in handling 
and marking, a common experience in entomo- 
logical work. 

The selection of which animals are killed, and 
which are sampled (as well as other events of a 
stochastic nature), is controlled by а pseudo- 
random number generator which supplies numbers 
within the population bounds. The appropriate 
process is then applied to the animal occupying 
that position in the population array “РОР” (after 
a check that that animal is in fact, ‘alive’). The 
pseudo-random number gencrator has Бсеп eare- 
fully tested, and supplies a random sequence of 
real numbers, uniformly distributed in the long 
run between 0 and 1. This number is multiplied 
by the standing population size at the time to 
obtain the order number of the animal selected. 
The experimenter supplies a set of six random 
numbcrs, obtaincd from one of the published 
tables (e.g. Rohlf & Sokal 1969) to seed the 
generator separately for each process. 

In each population cycle the population under- 
goes a death phase first, followed by a birth phase. 
The numbers of dcaths and births are specified 
by separate sub-routines according to the experi- 
menter’s instructions. 

DEATH: In the death phase a random number 
is generated and this animal is examincd. If it 
has a value of 99 its death is recorded; if its 
value indicatcs it is a marked animal its mark is 
recorded in another array ‘WLDKIL’; if its value 
is zero, this indicates the animal in this array 
position is dead, and it is ignored for the moment. 
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In the first two cases the value іп the array “РОР” 
is set to zero, and another random number is 
obtained. This procedure continues until the re- 
quired number of deaths has oecurred. 

BirTH: In the subsequent birth phase the array 
‘POP’ is checked through and zero spaces are 
changed to 99 (new-born animals) until the speci- 
fied number of births has occurred. If there are 
insufficient zero spaces in array “РОР”, its upper 
bound is inereased with new-born animals. 
SAMPLING: The array “РОР” is then sampled. 
Animals are again selected by the random number 
generator; each is examined and, if dead, another 
is seleeted. If the animal has not been previously 
captured (i.e. equals 99) its value is set to 1. 
Previously marked animals thus have marks other 
than О or 99—at the start of each cycle, all pre- 
viously marked animals have their value multiplied 
by 10 so an animal caught in the previous cycle 
would, on capture, have a mark of 10. A 1 is 
then added, and this animal's mark is changed to 
11. Similarly, an animal earrying a mark 10101 
would have been marked in the present, last but 
one and last but three cycles. The marks of such 
previously captured animals (recaptured during 
sampling) are noted іп an array ‘MARX’, for 
later printout, and the mark is analysed by the 
main program (for the Lincoln Index and Зећпа- 
bel estimators) and by the sub-routines BAILEY 
and JOLLY for the Bailey Triple Catch estimator 
(Bailey 1951, 1952) and Jolly's Stochastie esti- 
mator (Jolly 1965). 

Finally, a random number is drawn; if this is 
less than the specified probability of an animal 
being injured during marking, this animal's mark 
is transferred to another array ‘SMPLDD’, and its 
position in the population array “РОР” is set to 
zero (animal dead). This cycle is repeated until 
the required size of sample has been obtained. 

This sequence of deaths, births and sampling 
is thus carried out on the population for as many 
time periods as the experimenter may select. Such 
a sequenee, representing a population over a 
period of time, is herein referred to as an experi- 
ment. For each set of paramaters—seed value, 
birthrate, deathrate, sample size, number of time 
periods—the experimenter may also nominate 
the number of times he wishes this experiment 
repeated independently (or replicated). This is 
done on the parameter card (see facsimile of print- 
out). For each rcplieate, the experimenter supplies 
a new set of six random numbers, as mentioned 
above, which seed the pseudo-random number 
generator for probabalistic death or birth fune- 
tions, random selection of animals for death and 
capture, variation in sample size, and whether an 
animal chances to be killed during marking. 


Finally, a sequence of different experiments may 
be earried out in a single run. 
The program used here consists of a main pro- 
gram which acquires data (the parameters and 
random number sceds) and earries out the actual 
manipulation of the population. There аге a series 
of sub-routines and functions, as follow: 
Sub-routine SAMPLE—caleulates sample size 
BIRTH —caleulates number of 
new animals in each cycle 
DEATH —caleulates number of 
deaths in each cycle 
RNUM-—supplies random num- 
bers, n, (0 n x 1) 
LINC—estimates the population 
size by the Petersen or 
Lincoln Index method 


” 


33 
Function 


Sub-routine 


m SCHNAB—estimates the population 
size by Schnabel’s method 
BAILEY—estimates the population 


size Бу Bailey’s 'triple- 
catch’ method 
JOLLY—estimates the population 


size by Jolly’s stochastic 


method 
Ая LINCSM—collates the results of the 
successive Petersen esti- 


mates over the several re- 
plicates, supplying a mean 
and standard deviation for 


each cycle 

55 SNBSUM —collates results for 
Schnabel estimator as 
above 

T BALSUM —collates results for Bai- 


ley cstimator as above 

Я JOLSUM —collates results for Jolly 
estimator as above 

5 PRINT—organizes output from the 
four collating sub-routines 

РА TRUVAL—prints ош summary of 
true values of populations 
at the end of each sequence 
of experiments 

A faesimile of the main program and the sub- 
routines (excluding the last six, book-keeping 
sub-routines) is reproduecd as an Appendix to 
this paper. 

Output is produced on the line printer for each 
experiment. After identification, all the parameters 
for the experiment are given. Optionally, for each 
cycle, there is a summary of population size, num- 
ber of births and deaths, and size of sample 
drawn, and the marks of marked animals caught 
in samples, of those killed in marking, and of 
those killed ‘in the wild’ are displayed. (As a 
further option, the whole population may be 
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printed out also, but this obviously consumes a 
lot of paper and is only used in special cases.) 
The number of animals marked and unmarked, 
and the percentage marked, is printed, for each 
cycle. 

At the end of each experiment, the four sets of 
estimates of the population, together with other 
relevant statistics (e.g. cstimated standard devia- 
tion) is printed. The unbiased forms of the esti- 
mators are used to obtain population size (Cor- 
mack 1968). 

After all replicates of the experiment have bcen 
completed, summaries of the estimates for cach 
model are printed, and a table of the actual popu- 
Jation values in each replicate. 


RESULTS 

Thirteen separate expcriments were performed; 
each had different birth and death rates, as 
summarized in Table 1. Each experiment was 
replicated 20 times. The starting value for each 
population was 1,000, and the sample size 100 + 
10. The probability that an animal was injured 
during capture and marking was set at three pcr 
cent. Thus experiment 1 to 3 were with stable 
populations (except for the slight loss duc to in- 
jury) where birth and death rates were equal. In 
experiments 4 to 7 the populations undergo simple 
growth or decay; in 8, 10 and 12 the populations 
go through a growth and then a decay cycle; and 
in 9, 11 and 13, the sequence is reversed. 

The results for each expcriment are summarized 
graphically (Figs. 1-13). The upper half of each 
graph shows the results obtained in two replicates, 
chosen at random, by each of the four estimator 
models. Below each of these the mean + one 


TABLE 1 
Experiment! Periods Birthrate Deathrate 
1 All 0 0 
2 АП 0.10 0.10 
3 АП 0.20 0.20 
4 АН 0.10 0.0 
5 АП 0.0 0.10 
6 АП 0.20 0.10 
7 АП 0:10 0.20 
8 1-5 0-10 0.0 
6-10 0.0 0.10 
2 1-5 0.0 0.10 
6-10 0.10 0.0 
10 1-5 0.20 0.0 
6-10 0.0 0.20 
11 1-5 0.0 0-20 
6-10 0.20 0.0 
12. 1-5 0-20 0.10 
6-10 0.10 0.20 
1З 1- 0.10 0.20 
6-10 0.20 0.10 


standard deviation for all twenty replicates is 
shown for each sampling period. The use of the 
standard deviation as a measure of spread is 
justified if all replicates are plotted (for example, 
as in Fig. 14, for Jolly estimator in experiment 1) 
when it will be seen that 60-70 per cent of points 
lie within the limits of = one standard deviation. 
Superimposed on each graph is the actual trend 
line of the population, for comparison. 

Fig. 15 shows the proportion of the population 
bearing marks at the end of each sampling period, 
for each experiment. 


DISCUSSION 


Of the four estimation modcls used, two 
(Lincoln Index and Schnabel) are based on the 
assumption that there are no gains (through birth 
and immigration) or losses (by death or emigra- 
tion) to the population, while the othcr two 
modcls purport to account for thesc. The mathe- 
matical background (including other assumptions) 
to these is discussed by Cormack (1968) апа 
Southwood (1966); both authors cite a large 
number of references to field application of 
these methods, the former being biased toward 
fisheries applications and the latter toward ento- 
mological studies. (Readers should note that South- 
wood has a number of errors in his mathematical 
formulae—marring an otherwise exccllent dis- 
cussion.) 

Such theoretical modcls allow biologists to 
obtain estimates of population size and other 
population attributes from mark-recapture studies, 
together with estimates of the standard deviations. 
Hopefully, this allows the biologist to place con- 
fidence intervals on his estimate(s) of the popula- 
tion: that is, he states that hc has, say, 95 per 
cent confidence that the true size of the population 
lies within (approximately) two standard devia- 
tions of the estimated value. Thus Cormack 
(1968) citcs examples of the use of Jolly’s (1965) 
stochastic method, and quotes percentage ‘accura- 
cies’ achieved. 

Models, however, are predicated on assump- 
tions, as are the statistical methods used to 
elucidate the models. Is it possible that the 
assumptions аге not justified, in which case cither 
the estimated population size, and/or the associa- 
ted variance estimate, could be misleading? To 
date, thc only way of testing this proposition was 
by the comparison of different models. If the 
answers did not correspond, the particular worker 
made a subjective judgment in favour of one or 
more models (e.g. Parr 1965). 

Manly (1970) approached this problem by 
computer simulation of a mark-recapture system. 
Unfortunately the method he used is not clearly 
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described іп the paper. Each population starts 
at zero and grows over a period of ten ‘days’. At 
the start of each day, births occurred (either 25, 
50, 75 or 100), and at the end of each day, deaths 
occurred, in accordance with one of twelve life 
tables. Sample size varied from approximately 10 
to 95 per cent. Thcre were 96 populations and 
each was sampled, independently, tcn times during 
each of the last five days giving ten independent 
sets of estimates of the population. Estimation 
models used were Fisher and Ford (1947), Jolly 
(1965) and Manly and Parr (1968). 

Manly analyses the results obtained in terms 
of ‘percentage rclative bias’ (і.е. [cstimate-truc 
value]/true value, as а pcrcentage) and пісап 
square error. He found that when its assumptions 
were met, the Fisher and Ford method gave re- 
sults similar to the Jolly method and would be 
useful whcn sample sizes wcre small; however, it 
assumes a constant birthrate which may be un- 
realistic in a field situation. He was concerned 
that Jolly’s method appeared to give standard 
error estimates correlated with the estimates to 
which they relate. Manly and Parr's method gave 
good rcsults when sample sizes were ‘moderately 
large’. Manly’s approach is interesting and inno- 
vative; regrettably, it docs not yicld any sort of 
Statistical tool that can be applied to other authors’ 
work. 

A different approach was uscd in the present 
study. Conccptually, the population ‘lives’ in an 
array of numbered sites in the computer. If an 
animal dies, its site becomes available to another, 
or new sites are created if песйс to house ап 
increasing population. By kecping account of how 
many sites are in usc (even if some are unoccu- 
pied), the program (by use of a pseudo-random 
number generator) ensurcs that animals are killed 
at random, and that all live animals have an equal 
chance of being sampled. Captured animals carry 
a mark that can be translated to yield their com- 
plete trapping history. A facility to allow for thc 
possibility that a sampled animai is killed in 
handling was included, firstly bccause this is realis- 
tic, and secondly becausc Jolly's modcl allows for 
this contingency. Sample size is not fixed, but 
varies between limits, which is analagous to the 
field situation where the size of sample is often 
dependent on mctcorological, human and logistical 
factors, opcrating around the biometrician's dcsign 
of the sampling program. 

Other options werc built into the program, so 
that, given a set of ficld parameters (approximate 
population size, birth and death rates, etc.), the 
population eould be simulated and various sam- 
pling regimes tried out. In addition, the four sub- 
routines used in the program to estimate popula- 


tions ean be detached and applied to actual 
capture-recapture data. 

An element of artificiality is produced by 
carrying through the entire death phase followed 
by the birth phase. As a result, in growing popula- 
tions some of the new-born animals are clustercd 
at the upper end of the array. А process to 
re-randomize the animals іп the array prior to 
sampling was tricd, but as this produced no 
difference in the population estimates (but appreci- 
ably increased the amount of computer timc used) 
this was not incorporated in the final modcl. This 
may not, in any event, bc important: Jolly (1965) 
points out that the important point is that the 
probability of capturing a marked animal must be 
the same as the probability of capturing any 
member of the population, and this is not changed 
by the arrangement used. 

The eventual aim of this project is to examine 
different population and sample size combinations, 
with scveral birth ratc/death rate combinations, to 
formulate sets of tables to yield more realistic 
confidence intervals for the models sclected. 


POPULATION ESTIMATES 


Figs. 1-13, read in conjunction with Table 1, 
are self-explanatory, and it does not seem neces- 
sary to discuss them in detail, Some general 
conclusions do emerge that are worth amplifying. 

1. Bailey's triple catch method never seems to 
yield particularly satisfactory results. In general, 
the estimates are not as good as those of the 
Lincoln index, and this model will bc disrcgarded 
in the discussion below. 

2. Schnabel's method yields excellent results in 
populations subject to neither birth nor death, or 
increasing (experiments 1, 4, 6 and, to a lesser 
extent, the growth phases of other experiments); 
it will be noted that thc standard deviation of 
the twenty replicates is uniform and reasonably 
compact (although the standard deviations csti- 
mated by the Schnabel model arc ludicrously 
small). This is a most unexpected conclusion as 
it has always been assumed that recruitment (as 
well as mortality) would be a source of error (e.g. 
Ricker 1958, Cormack 1968). Certainly the model, 
apparently, does not detect mortality: as will be 
seen in expcriments 2 and 3, the model predicts a 
growing population, whilc in 5 and 7 it predicts 
a more or less stable onc. Thus if a reliable (pre- 
ferably independent) cstimatc of birth and death 
rates is available, thc Schnabel model could pro- 
vide усгу good estimatcs of population size. 

3. Lincoln Index and Jolly's estimator generally 
perform best where the proportion of marked 
individuals is high (e.g. 5 and 7), in contrast to 
Schnabel, which performed best where the pro- 
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portion of marked individuals is low (Fig. 15); 
an apparent exception is seen in Fig. 1 where 
prediction by Schnahel was good and the per- 
centage marked was high, but with neither hirth 
nor death the Schnabel model becomes an ex- 
tended Lincoln estimate. Ricker (1945) claims 
that the Schnabel estimator has maximum effi- 
ciency when a negligible proportion is marked. 
In this context, note that the Schnabcl estimator 
did not perform as well in experiment 4 as in 
experiment 6, where the percentage marked was 
lower (Fig. 15). This is also illustrated in the 
experiments where the population is declining. For 
the present regime, where sample sizc is about 10 
per cent or less of the population, thcre appears 
to be little to choose hetween Lincoln Index and 
Jolly’s estimator. However, Jolly's performance is 
markedly hetter whcre birth and death rates are 
higher—the succession of Lincoln Index estimates 
show no improvement in performance with time, 
hut the standard deviation for the Jolly cstimates 
is smaller in the last pcriods than in the first few. 
(This is most clearly seen in expcriment 3.) In 
addition, the Jolly model provides estimates of 
hirth and death rates (or, more broadly, dilution 
and loss rates (Southwood 1966) (۰ 

4. The variance (or standard deviation) esti- 
mates given by the Jolly model are, as Manly 
(1970) points out, loosely correlated. with the 
estimate to which they apply. This was observed 
in the present study also, but it is not really sur- 
prising. When the true population size is known 
it is observed that the standard deviation asso- 
ciated with cxcessively large estimates usually 
encloses the true population value. In fact, the 
standard deviations are highly realistic from this 
point of view, although this would be little con- 
solation to a biologist endeavouring to estimate 
the population size within a run of several such 
large estimates of size and standard deviation. 

5. It will be observed that thc mean of the 
twenty replicates for all the models follows (with 
certain exceptions in Schnabcl) the true popula- 
tion curve rather wcll. This suggests that im- 
provements in estimates could he obtained hy 
sub-sampling from populations, where possible, 
and averaging the estimates rather than relying 
on point estimates. 

Further analyses of the present data are as yet 
incomplete, for example the relationship hctween 
estimates and standard deviations, loss and dilution 
rates derived by the cstimation models. The prc- 
sent study will he broadened hy examining the 
models with different initial population sizes, by 
using different sampling rates, and hy allowing 
greater variation in sample size, with constant 
birth and death rates as in experiments 1 to 7. 


These should lead to further generalization con- 
cerning the validity of the various models. 


CONCLUSIONS 


For a series of estimates in a mark-recapture 
study, Jolly's (1965) model gives the most satis- 
factory results in situations similar to those tested 
here, as it includes realistic standard deviation 
estimates and other data. For the special case of 
a stable closed population, or one undergoing 
growth, the Schnabel estimator gives excellent 
results, although its standard deviation estimate is 
totally misleading. For those cases where multiple 
or individual marking of animals is difficult (as 
in many insects) a serics of Lincoln Index esti- 
mates gives quite satisfactory results. 
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DESCRIPTION OF FIGURES 1-16 


Fics. 1-13—-These figures summarize the results of experiments 1-13 respectively (see life- 
history parameters in Table 1). In all cases, the starting value for the population was 1000; 
sample size was 100 + 10; and the probability of loss during handling was 3%. 

Thc upper half of each figure shows the actual estimates obtained in periods 2-9 in two 
replicates selected randomly, using (left to right) the Lincoln, Schnabel, Bailey and Jolly 
models. Where estimates excecd the limits of the scale (shown on the right) the estimate is 


shown in numbers. 


The lower half of each figure (ІсЁ: hand scale) shows the mean (horizontal bar) = 1 
standard deviation (vertical bar) over the twenty replicates, using the same тойс]. 
The solid line shows the actual population trend. 


Fic. 14—Distribution of estimates by Jolly model for births = deaths = 0. Each point 
represents an estimate in the appropriate pcriod, the horizontal bar the mean, and the vertical 


bar + 1 standard deviation. 


Fic. 15—Percentage of marked individuals, by timc periods, for each of the experiments 1-13. 


Ею. 16—Flow diagram of main program. 
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